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PREFACE 


ORACLE  (Oversight  of  Resources  and  Capability  for  Logistics 
Effectiveness)  is  a  methodology  that  was  developed  to  abstract  aggregate 
relations  from  the  U.S.  Air  Force's  Recoverable  Item  Requirements 
Computation  System,  which  is  also  known  as  D041.  Because  the  current 
method  for  computing  the  ORACLE  database  is  designed  to  work  with  only  a 
simplified  version  of  D041,  the  purpose  of  Volume  1  of  this  Note  is  to 
investigate  how  this  method  can  be  modified  to  work  with  the  actual 
version.  This  volume  presents  technical  details,  and  it  is  intended  for 
readers  who  may  wish  to  implement  or  extend  this  methodology,  or  to 
understand  its  technical  limitations.  Volume  II  examines  how  an 
ORACLE-like  procedure  could  be  used  to  improve  requirements  forecasts. 

This  Note  was  prepared  for  the  Office  of  the  Assistant  Secretary  of 
Defense  for  Production  and  Logistics  under  the  auspices  of  RAND's 
National  Defense  Research  Institute,  the  OSD-sponsored  Federally  Funded 
Research  and  Development  Center  at  RAND.  The  work  was  performed  under 
the  project  titled  "Effective  Modeling." 
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SUMMARY 


.'ORACLE  (Oversight  of  Resources  And  Capability  for  Logistics 
Effectiveness)  is  a  methodology  that  assesses  the  effects  of  changing 
certain  resource  levels  on  the  peacetime  materiel  readiness  and  wartime 
sustainability  of  U.S.  air  forces,  to  improve  the  estimation  and 
justification  of  resource  requirements.  The  Air  Force  Logistics  Command 
(AFLC)  uses  the  D041  system  to  help  manage  approximately  150,000 
aircraft  components.  During  each  quarter,  the  D041  system  is  used  to 
estimate  how  much  of  each  component  should  be  repaired  and  purchased  for 
about  three  years  into  the  future.  The  goal  of  the  ORACLE  methodology 
is  to  construct  an  aggregate  database  having  the  following  features: 

The  database  is  an  additional  product  of  a  standard  D041  quarterly 
computation,  it  is  small  enough  to  fit  in  a  portable  microcomputer,  and 
it  can  be  readily  manipulated  by  a  spreadsheet-like  program  to  mimic  in 
aggregate  form  the  responses  of  D041  to  program  changes. 

The  current  ORACLE  methodology  was  designed  to  work  with  only  a 
simplifed  version  of  D041,  not  the  actual  version.  D041  is  now  being 
modified  to  incorporate  the  logic  that  is  in  the  Logistics  Management 
Institute 1 s f(LMf)  Aircraft  Availability  Model  (AAM) .  In  particular, 

D041  will  incorporate  the  AAM1  s  approach  to  common  items,  indenture 
relationships,  and  an  availability  objective  based  upon  no 
cannibalization.  Because  the  current  methodology  does  not  allow  for  any 
of  these  new  features,  this  Note  investigates  how  these  features  can  be 
included. 

Altogether,  five  different  methods  for  developing  the  ORACLE 
database  are  described.  The  first  method  is  based  upon  differentiating 
the  necessary  conditions  for  optimality  and  mean  pipeline  equations, 
which  is  an  extension  of  the  approach  used  by  the  preliminary  version  of 
ORACLE.  Unfortunately,  the  indenture  relationships  make  this  approach 
difficult  to  implement.  The  second  method  is  an  approximate  procedure 
in  which  the  indenture  relationships  are  essentially  eliminated.  The 
third  method  retains  the  indenture  relationships  but  in  a  simplified 


form.  The  fourth  is  based  upon  making  multiple  runs  of  the  AAM  for  the 
same  year.  The  fifth  method  uses  regression  analysis  to  analyze  data 
obtained  for  several  years.  These  methods  are  compared  with  respect  to 
simplicity,  accuracy,  and  execution  time. 

The  only  method  for  which  there  is  some  computational  experience  is 
the  final  one,  regression  over  time.  Volume  II  discusses  the 
performance  of  that  approach  in  the  context  of  requirements  forecasting. 
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I.  INTRODUCTION 


ORACLE  (.Oversight  of  Resources  and  Capability  for  Logistics 
Effectiveness)  is  a  methodology  that  was  developed  to  relate  dollars 
expended  on  recoverable  components  to  the  goals  set  in  the  Planning, 
Programming,  and  Budgeting  (PPB)  process.  [1]  This  methodology  is 
designed  to  systematically  abstract  aggregate  relations  between  dollars 
and  goals  from  the  U.S.  Air  Force's  Recoverable  Item  Requirements 
Computation  System,  which  is  also  referred  to  as  D041.  However,  the 
current  version  of  ORACLE  is  designed  to  work  with  only  a  simplified 
version  of  D041,  rather  than  the  actual  version. 

The  Program  Objective  Memorandum  (POM)  describes  the  activities  and 
capabilities  to  be  achieved  over  a  five-year  period,  as  well  as  the  time- 
phased  resources  required  to  achieve  them.  The  Air  Force  Logistics 
Command  (AFLC)  uses  the  D041  computation  as  the  basis  for  its  inputs  to 
the  POM.  However,  the  D041  computation  makes  forecasts  for  only  about 
three  years  beyond  the  asset  cutoff  date,  while  the  POM  years  extend  to 
about  seven  years  beyond  asset  cutoff.  In  the  past,  a  cost-per-f lying- 
nour  rate  has  been  used  to  forecast  POM  requirements  beyond  D04l's  horizon 
but  in  recent  years  these  forecasts  have  been  substantial  underestimations 

This  Note  will  investigate  the  following  two  questions: 

1.  How  can  the  ORACLE  estimation  methodology  be  modified  so  that 
it  can  work  with  the  actual  version  of  D041? 

2.  How  can  an  ORACLE- like  procedure  be  used  to  improve  the 
accuracy  of  the  POM  requirements  forecasts? 

Volume  I  will  consider  the  first  question,  and  Vol.  II  will  consider  the 
second . 

To  develop  the  ORACLE  dst-obzst;  for  a  weapon  system  (or  end  item) 
included  within  D041,  the  first  step  is  to  identify  a  set  of  independent 
variables  Z^,  m  =  1,  ...,  M  (such  as  an  availability  target  and  total 
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flying  hours)  and  a  set  of  dependent  variables  D  ,  n  =  1,  ....  S'  (such 

n 


as  buy  requirement,  serviceable  assets,  and  base  repairs).  The  second 
step  is  to  run  D041  for  a  nominal  case,  which  means  for  a  specified  set 
of  values  for  the  independent  variables  for  each  weapon  system.  And  the 
final  step  is  to  compute  the  partial  derivative  of  each  dependent 
variable  with  respect  to  each  independent  variable,  as  evaluated  at  the 
nominal  case.  The  ORACLE  database  for  each  weapon  system  consists  of 


the  collection  of  partial  derivatives  3D  /3Z  ,  where  m  =  1,  ...,  M  and 

n  m 


n  =  1 ,  .  .  .  ,  N . 

After  the  ORACLE  database  has  been  constructed,  the  outputs  of 


D041  can  be  approximated  in  the  following  way.  For  a  given  weapon 
system,  let  Z  ^  be  the  value  of  the  mth  independent  variable  for  the 


0 


nominal  case,  and  let  D  be  the  corresponding  value  of  the  nth 


dependent  variable  as  computed  by  D041.  If  a  scenario  is  considered  for 


which  the  independent  variables  are  Z^,  for  m  =  1,  ...,  M,  then  the 


corresponding  value  of  the  nth  dependent  variable  that  would  be  computed 
by  D041  can  be  approximated  with 

M 

0 

+  7  oL)  /d  L  it 

m  m 

rti=l 


D  =  D  °  +  y'  3D  / 3Z  (Z  -  z  °)  . 
n  n  £-J  n  mm  m 


(1) 


D041  is  now  being  modified  to  incorporate  the  logic  that  is 
currently  in  the  Logistics  Management  Institute's  (LMI)  Aircraft 
Availability  Model  (AAM).[2]  In  particular,  D041  will  incorporate  the 
AAM's  approach  to  common  items,  indenture  relationships,  and  an 
availability  objective  based  upon  no  cannibalization.  Because  the 
current  method  for  computing  the  ORACLE  database  does  not  allow  for  any 
of  these  features,  one  purpose  of  this  Note  is  to  investigate  how 
these  new  features  could  be  included.  Next,  each  of  these  features  is 
described  separately. 

The  term  "common  item"  refers  to  an  item  that  is  installed  on  more 
than  one  weapon  system.  Thus  when  a  common  item  is  backordered,  the 
availabilities  of  several  weapon  systems  may  be  degraded.  When  a  spare 
unit  of  a  common  item  is  procured,  however,  several  weapon  systems  share 
the  benefit.  Based  upon  the  relative  amount  of  flying  hours,  the  approach 


taken  by  the  AAM  is  to  allocate  a  portion  of  a  common  item  s  backorders 
and  cost  to  each  associated  weapon  system. 

The  indenture  relationships  indicate  which  items  are  subassemblies  of 
other  items.  A  Line  Replaceable  Unit  (LRU)  refers  to  a  component  that  can 
be  removed  and  replaced  on  the  flight  line;  a  Shop  Replaceable  Unit  (SRU) 
refers  to  a  subassembly  that  can  be  removed  and  replaced  in  a  repair  shop. 
The  AAM  treats  the  effect  of  LRU  backorders  differently  from  that  of  SRU 
backorders.  Although  an  LRU  backorder  affects  aircraft  availability,  an 
SRU  backorder  lengthens  LRU  repair  times,  thus  lowering  the  probability  of 
an  LRU  spare  being  in  a  serviceable  condition. 

The  objective  function  used  by  the  AAM  is  the  probability  that  a 
random  aircraft  is  not  missing  any  of  its  reparable  LRUs,  assuming  that 
there  is  no  cannibalization.  In  particular  it  is  assumed  that  the  holes 
in  the  aircraft,  because  of  a  lack  of  serviceable  LRUs,  are  uniformly 
spread  over  all  appropriate  aircraft.  This  objective  function  depends 
upon  the  expected  numbers  of  LRU  backorders  at  the  base  echelon,  which  in 
turn  depend  upon  the  expected  numbers  of  LRU  backorders  at  the  depot 
echelon  and  upon  the  expected  numbers  of  SRU  backorders  at  both  the  base 
and  depot  echelons. 

Section  II  will  state  the  equations  in  the  AAM,  as  well  as  the 
underlying  assumptions  and  approximations.  This  model  has  the  structure 
of  a  constrained  nonlinear  programming  problem,  and  it  is  solved  using  a 
marginal  analysis  technique. 

Section  III  will  present  a  method  for  computing  the  ORACLE  database 
for  the  AAM  based  upon  differentiating  the  necessary  conditions  for 
optimality  and  mean  pipeline  equations,  which  is  an  extension  of  the 
approach  used  by  the  current  version  of  ORACLE.  Because  only  LRUs  are 
considered  in  the  current  version,  that  version  is  able  to  compute  the 
ORACLE  database  by  treating  only  one  item  at  a  time.  However,  one  effect 
of  including  indenture  relationships  is  to  couple  the  partial  derivatives 
for  an  LRU  with  those  for  all  associated  SRUs .  In  that  case,  computing 
the  ORACLE  database  requires  solving  a  system  of  simultaneous  linear 
equations  for  each  LRU  in  the  weapon  system. 


mmm, 


Because  the  indenture  relationships  make  the  foregoing  approach 
difficult  to  implement,  Sec.  IV  will  discuss  four  alternative 
estimation  procedures:  The  first  method  treats  the  SRUs  as  though  they 
were  LRUs,  which  essentially  eliminates  the  indenture  relationships;  the 
second  retains  the  indenture  relationships  but  in  a  simplified  form;  the 
third  is  based  upon  making  multiple  runs  of  the  AAM  for  the  same  year; 
and  the  fourth  uses  regression  analysis  to  analyze  data  obtained  for 
several  years. 

Altogether,  Secs.  Ill  and  IV  will  present  five  methods  for 
constructing  the  ORACLE  database.  Section  V  will  compare  these  methods 
with  respect  to  simplicity,  accuracy,  and  execution  time. 


II.  AIRCRAFT  AVAILABILITY  MODEL 


This  section  will  describe  the  formulas  and  assumptions  underlying 
the  Aircraft  Availability  Model  (AAM)  based  upon  the  most  recent 
documentation  for  that  model. [2]  Because  many  of  these  formulas  were 
not  explicitly  stated  in  this  documentation,  it  was  necessary  to  verify 
their  accuracy  by  checking  with  the  developers  of  the  AAM. 

Consequently,  this  section  may  provide  the  most  complete  mathematical 
statement  of  the  AAM  that  is  available. 

The  following  simplifying  assumptions  will  be  used: 

1.  There  is  only  one  level  of  indenture,  which  means  that  each 

SRU  is  a  subassembly  of  an  LRU,  rather  than  of  another  SRU; 

2.  An  SRU  can  be  a  subassembly  for  only  one  type  of  LRU. 

Although  AAM  does  not  require  either  of  these  assumptions,  they  will  be 
used  to  simplify  the  presentation  given  here. 

OBJECTIVE  FUNCTION 

The  objective  function  of  AAM  will  be  written  as  a  direct  function 
of  the  expected  numbers  of  LRU  backorders  at  bases.  However,  it  will  be 
shown  later  that  the  expected  numbers  of  LRU  backorders  at  bases  depend 
upon  the  expected  numbers  of  LRU  backorders  at  the  depot  echelon  and 

upon  the  expected  numbers  of  SRU  backorders  at  both  the  base  and  depot 

echelons . 

A  weapon  system  refers  to  a  particular  mission  and  design  of  an 
aircraft,  which  may  have  more  than  one  series.  For  a  given  weapon 
system,  the  index  i  will  refer  to  the  ith  LRU  and  the  index  k  will  refer 
to  the  kth  series.  The  index  e  will  refer  to  the  repair  echelon:  e  =  1 
represents  the  base  echelon  and  e  =  2  the  depot  echelon. 

Suppose  that  a  particular  weapon  system  has  been  specified  and 
define 

akiO  =  (luantity  Per  application  of  LRU  i  on  series  k 


b,  .  =  fraction  of  aircraft  in  series  k  that  includes  LRU  i 
ki 


n.  =  number  of  aircraft  in  series  k. 
k 

Thus  the  number  of  units  of  LRU  i  installed  on  aircraft  of  series  k  is 
tkiO  aki0^kink’ 

If  K  is  defined  to  be  the  number  of  series  in  the  specified  weapon  system, 
then  the  number  of  aircraft  in  that  weapon  system  is 
K 

N  =  £  v 

k=l 

Define 

=  expected  number  of  backorders  for  LRU  i  at  echelon  e. 

Thus  EBCh  refers  to  the  expected  number  of  backorders  for  LRU  i  at  all 
bases.  Because  the  AAM  allows  the  possibility  that  LRU  i  could  be 
installed  upon  several  series  of  the  same  weapon  systems,  or  upon 
different  weapon  systems,  it  is  necessary  to  compute  the  fraction  of  the 
total  backorders  that  is  associated  with  a  particular  series  of  a 
particular  weapon  system,  which  is  done  as  follows.  First  define 

=  number  of  flying  hours  for  series  k  of  the  specified  weapon 
system 

IP^  =  number  of  flying  hours  accumulated  by  units  of  LRU  i  over 
all  weapon  systems 


and  then  compute  the  use  factor 


7 


i 


■§ 

m 


The  AAM  considers  U^EIIO  ^  to  t^le  exPectec^  number  of  backorders  for 
LRU  i  that  is  associated  with  series  k  of  the  specified  weapon  system. 

Assuming  that,  the  backorders  are  randomly  distributed  over  the 
available  aircraft,  then  the  probability  that  an  aircraft  of  series  k  is 
not  missing  a  unit  of  LRU  i  is 


q.  .  =  1  -  b.  .  +  b.  .(1  -  U.  .EBO.ni/t.  ,A)akiO. 

kiO  ki  ki  ki  iOI  kiO 


The  probability  that  a  randomly  chosen  aircraft  in  the  specified 
weapon  system  is  not  missing  a  unit  of  LRU  i  is 


qio  ■  £  <vN)w 

k=l 

If  I  is  defined  to  be  the  number  of  LRUs  in  the  weapon  system  and 
it  is  assumed  that  the  various  LRUs  fail  independently,  then  the 
probability  that  a  randomly  chosen  aircraft  is  not  missing  any  of  its 
reparable  LRUs  is 


*  ■  n  «io- 


(2) 


i=l 

This  expression  is  a  measure  of  aircraft  availability  without 
cannibalization  and  is  the  objective  function  that  will  eventually  be 
maximized . 

EXPECTED  SRU  BACKORDERS  AT  DEPOT 

The  index  pair  (i,  j)  will  be  used  to  refer  to  the  jth  SRU 
belonging  to  the  ith  LRU.  However,  if  j  =  0,  then  (i,  0)  will  refer  to 
the  ith  LRU.  Define 


X..  =  daily  demand  rate  for  SRU  (i,  i)  at  echelon  e 

lje  J 


D(Tj  =  depot  condemnation  rate  for  SRU  (i,  j) 


DL„  =  depot  procurement  time  for  SRU  (i,  j) 


8 


RT . .  =  average  repair  time  for  SRU  (i,  j)  at  echelon  e 

t  J  e 

x..  =  stock  level  for  SRU  (i,  j)  for  each  site  at  echelon  e 

ije  J 

m. .  =  mean  pipeline  for  SRU  (i,  j)  for  each  site  at  echelon  e. 

ije  J 

The  depot  pipeline  for  an  SRU  consists  of  the  units  being  repaired 
at  the  depot,  plus  those  on  order  from  a  commercial  vendor.  Thus  the 
mean  depot  pipeline  for  SRU  (i,  j)  is 


m..0  =  X .  [ RT . ( 1  -  DC..)  +  DL..DC..], 
ij2  ij  2 1  ij2  ij'  ij 


(3) 


where  e  =  2  refers  to  the  depot  echelon. 
Define 


p(y|m..  )  =  probability  of  having  y  units  in  the  pipeline  given 
ije 

that  the  mean  is  m. ,  . 

ije 

The  AAM  uses  either  a  Poisson  or  negative  binomial  distribution  as 
the  probability  distribution  for  each  pipeline.  Given  the  stock  level 
Xij2’  t^e  exPectec*  backorders  for  SRU  (i,  j)  at  the  depot  are  given  by 

EB0--  =  ^  (y  -  x. .7)p(y|m. .,) . 


ij2 


ij2 


ij2J 


y>Xij2 


EXFECTED  LRU  BACKORDERS  AT  DEPOT 

The  depot  pipeline  for  an  LRU  consists  of  the  units  on  order,  plus 
those  being  repaired  at  the  depot,  plus  those  awaiting  SRUs  at  the  depot 
before  being  repaired.  However,  AAM  makes  the  approximation  that  the 
calculation  of  LRU  depot  backorders  can  be  uncoupled  from  that  for  SRU 
depot  backorders.  In  particular,  the  mean  depot  pipeline  for  LRU  is 
computed  as 


m . 


i02 


=  Xi02[RTi02(1  '  DC^  +  DL^DC,J, 


iO' 


iO  iO 


(4) 


which  is  the  same  as  Eq.  (3)  except  that  j  =  0. 

Later  in  Eq.  (6),  the  mean  SRU  backorders  at  the  depot,  which  should 
be  affecting  LRU  backorders  at  the  depot,  are  added  to  the  mean  LRU 
pipeline  at  the  base  echelon.  In  effect,  AAM  assumes  that  each  additional 
SRU  depot  backorder  will  cause  an  additional  LRU  to  be  awaiting  parts  at 
the  depot,  causing  an  additional  LRU  to  be  backordered  at  the  depot, 
adding  one  unit  to  the  mean  LRU  pipeline  at  the  base  echelon.  This 
approximation  will  overestimate  the  mean  LRU  base  pipeline  for  two 
reasons:  A  single  LRU  may  be  awaiting  parts  for  more  than  one  SRU,  and  an 
additional  LRU  that  is  awaiting  parts  may  not  necessarily  cause  an 
additional  LRU  to  be  backordered. 

Given  the  stock  level  x^^  an<^  ignoring  the  SRU  contribution,  the 
expected  backorders  for  LRU  i  at  the  depot  is 


EBOi02=Z  (y  ’  xi02)p(yl”i02)' 


y>x 


i02 


EXPECTED  SRU  BACKORDERS  AT  BASES 

Define 


RTS„  =  fraction  reparable  this  station  (at  base  level) 
for  SRU  (i,  j) 


OST^j  =  order  and  ship  time  from  depot  to  base  for 


SRU  (i,  j) 


OIM^j  =  fraction  of  depot  demand  for  SRU  (i,  j) 
originating  from  bases 


NB .  =  number  of  activities  that  use  LRU  i  at  echelon  e. 


The  AAM  treats  the  depot  echelon  as  though  it  is  a  single 
aggregated  entity,  which  corresponds  to  having  NB^  =  1.  However, 
multiple  bases  are  allowed,  which  corresponds  to  having  NB^  >  1. 

The  base  pipeline  for  an  SRU  consists  of  the  units  on  order  from 
the  depot,  plus  those  being  repaired  at  the  base.  However,  the  units 
on  order  must  include  the  requisitions  that  have  been  backordered  at  the 
depot,  as  well  as  those  in  transit  between  the  depot  and  base.  Because 
OIM^  is  the  fraction  of  depot  backorders  that  correspond  to  orders 
received  from  the  base  echelon,  the  mean  pipeline  for  SRU  (i,  j)  at  an 
average  base  is  computed  as 


m...  =  { X  RTS . . RT  +  X..,(l  -  RTS..)OST.. 
ijl  ijl  ij  ijl  ijl  ij  ij 


+  OIM . .EBO.  .  „ } /NB  .  , 
ij  ij  2  il’ 


and  the  expected  backorders  at  all  bases  are 


Eaoiji  -  E 

y>xiji 


NBn(y  -  XjjppCyl.^p. 


EXPECTED  LRU  BACKORDERS  AT  BASES 

The  base  pipeline  for  an  LRU  consists  of  the  units  on  order,  plus 
those  being  repaired  at  the  base,  plus  those  awaiting  parts  at  the  base 
before  being  repaired.  The  units  on  order  include  the  requisitions 
that  have  been  backordered  at  the  depot,  plus  those  in  transit  between 
the  depot  and  base. 

The  AAM  makes  two  key  assumptions  when  computing  the  mean  LRU  base 
pipeline:  One  LRU  at  a  base  is  awaiting  parts  for  each  backordered  SRU 

at  that  base,  which  ignores  the  possibility  that  a  single  LRU  could  be 
waiting  for  several  SRUs ;  and  the  fraction  (1  -  OIM„)  of  the  backordered 
SRUs  at  the  depot  will  have  the  same  effect  as  the  same  number  of 
backordered  LRUs  at  the  depot,  which  enables  the  calculation  of  depot 
SRU  backorders  to  be  uncoupled  from  that  for  depot  LRU  backorders.  For 
simplicity,  this  section  has  also  assumed  that  each  SRU  can  be 
associated  with  only  one  type  of  LRU.  Under  these  circumstances,  the 
mean  pipeline  for  LRU  i  at  an  average  base  is 


"iOl  =  lXi01RTSi08Ti01  +  XiOl‘!  -  RTSiO)OSTiO 


+  EBOi02  +  £  IEBOiil  +  (1  -  °IVEBOij21)/NBil- 


where  J(i)  is  the  number  of  SRUs  associated  with  LRU  i.  Thus  the  expected 


backorders  for  LRU  i  at  all  bases  are 


>i01  -  £  NBil(y  •  xi01)p(y,mi01)- 


y>Xi01 


SAFETY  STOCK  CONSTRAINTS 


The  safety  stock  is  defined  to  be  the  difference  between  the  stock 


level  and  the  mean  pipeline;  that  is,  x..  -  m..  .  Although  not 

ije  lje 

mentioned  in  the  AAM  documentation,  the  planned  D041  implementation  of 


the  AAM  will  restrict  the  safety  stock  to  be  nonnegative.  For 
completeness,  this  section  will  allow  an  upper  bound  restriction  to  be 
present  as  well.  Define 


S(m..  )  =  pipeline  standard  deviation  for  SRU  (i,  j) 
i  J  e 


at  echelon  e  as  a  function  of  the  corresponding  mean. 


For  example,  D041  currently  assumes  that 


SO,..  )  =  a1/20»..  )<B+I)/2, 

ije  ije 


where  a  and  6  are  constants  determined  by  a  statistical  study.  For  the 


purposes  of  this  section,  it  is  assumed  that  the  stock  level  x..  has 

ije 

the  constraint 


LK. .  <  (x.  .  -  m. .  )/S(m . .  )  <  UK . . 

ije  ije  ije  ije  ije 


for  i  =  1,  ...,  I,  j  —0f  ....  J(i),  and  e  =  1,  2,  where  LK..  and  UK.,  a: 

ije  ije 

specified  lower  and  upper  bounds.  For  example,  LK^e  =  0  corresponds  to 
the  planned  norncgativity  restriction  in  D041. 
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OPTIMIZATION  PROBLEM 

The  final  task  is  to  formulate  the  optimization  problem.  First 
define 


STK .  .  =  available  assets  for  SRU(i,  j)  at  echelon  e. 

ije  J' 


The  available  assets  STK,.  include:  serviceable  spares  on  hand  and 

ije 


on  order,  plus  carcasses  waiting  for  induction  and  in  repair,  less  any 
backorders.  Because  broken  carcasses  are  included  as  available  assets, 
this  formulation  assumes  that  everything  that  can  be  repaired  will  be 
repaired.  However,  AAM  does  allow  another  option,  the  so-called  "repair 
option,"  in  which  the  magnitude  of  depot  repairs  is  also  a  decision 
variable . 

Next  define 


c„  =  replacement  cost  for  SRU(i,  j) 


Just  as  an  earlier  formula  prorated  backorders  for  common  items  among 
different  weapon  systems,  it  is  also  necessary  to  prorate  costs. 
Compute 


K  ) 

V.  J 
1 

akiObkiFk  ( 

k=l  ' 

which  is  the  proportion  of  the  total  item  flying  program  generated  by 


the  specified  weapon  system.  The  AAM  considers  V^c^  to  be  the  cost 
actually  charged  to  this  weapon  system  when  buying  one  unit  of 


SRU(i,  j). 


The  optimization  problem  is  to  determine  the  stock  levels  x..  for 

ije 

i=  1,  ...,  I,  j  =0,  ...,  J(i),  and  e  =  1,2  in  order  to  maximize  the 


objective  (2)  subject  to  the  safety  stock  restrictions  (7)  and  the 
budget  constraint 

I  J(i)  2 

ET  y  V  .  c  .  .  [NB  .  x..  -  STK..  ]+  <  B, 

l  ij  ie  ije  ijeJ 

i=l  j=0  e=l 


0 


otherwise 


v 


and  B  is  the  available  procurement  budget. 

The  AAM  has  the  structure  of  a  constrained  nonlinear  programming 
problem.  The  approximations  in  the  foregoing  equations  were 
deliberately  made  so  that  this  problem  can  be  efficiently  solved  by 
using  a  marginal  analysis  technique.  Because  the  depot  SRUs  are 
uncoupled  from  depot  LRUs,  it  is  possible  to  make  tradeoffs  among  SRUs 
for  a  given  LRU,  before  making  tradeoffs  among  LRUs.  The  LRUs  can  then 
be  ranked  in  terms  of  the  increase  in  the  objective  function  that  would 
occur  if  the  corresponding  spare  unit,  or  equivalent  investment  in  SRUs, 
were  added  to  the  inventory. 

For  a  component  that  is  common  to  several  weapon  systems,  the  AAM 
will  in  general  compute  a  different  stock  level  for  each  of  these  weapon 
systems.  In  this  circumstance,  the  AAM  will  implement  the  largest  of 
these  stock  levels,  so  that  the  weapon  system  "wanting"  the  component 
the  most  will  drive  the  decision.  For  the  purpose  of  computing  the 
ORACLE  database,  which  is  the  set  of  partial  derivatives  for  the 
dependent  variables  with  respect  to  independent  variables  belonging  to 
the  same  weapon  system,  this  interaction  between  different  weapon 
systems  will  be  ignored.  Incorporating  this  interaction  would 
correspond  to  adding  extra  terms  to  the  basic  linear  model  (Eq.  (1)),  where 
these  extra  terms  involve  the  partial  derivatives  for  the  dependent 
variables  belonging  to  one  weapon  system  with  respect  to  independent 
variables  belonging  to  other  weapon  systems.  Although  in  the  future  it 
may  prove  desirable  to  extend  ORACLE  so  that  these  extra  terms  are 
added,  such  an  extension  will  not  be  addressed  here. 
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III.  NECESSARY  CONDITIONS  FOR  OPTIMALITY 


For  the  purpose  of  developing  the  ORACLE  database  for  the  AAM,  it  is 
possible  to  use  formulas  given  by  Bigelow  [3]  for  computing  the 
derivatives  of  several  dependent  variables,  once  the  derivatives  for  the 
stock  levels,  base  repairs,  and  depot  repairs  are  known,  and  for  computing 
the  derivatives  for  base  and  depot  repairs.  The  only  formulas  missing  are 
for  computing  the  partial  derivatives  of  the  AAM  stock  levels  with  respect 
to  the  independent  variables.  This  section  will  show  how  to  compute  these 
latter  derivatives  by  differentiating  the  necessary  conditions  for 
optimality  and  mean  pipeline  equations,  which  is  an  extension  of  the 
method  originally  developed  by  B ige low . [ 3 , 4 ] 

A  key  approximation  made  throughout  this  section  is  to  use  the  normal 
distribution  as  the  probability  distribution  for  each  pipeline,  for  the 
formulas  for  computing  the  partial  derivatives  even  to  be  defined.  If 
discrete  distributions,  as  used  in  the  AAM  (the  Poisson  and  negative 
binomial)  were  used,  then  the  derivatives  would  not  be  defined.  It  is 
well  known  that  if  the  mean  pipeline  size  is  large,  then  the  normal 
distribution  does  provide  a  good  approximat ion . [ 5 ]  But  if  the  mean  is 
small,  the  normal  approximation  introduces  some  error.  However,  this 
approximation  is  not  needed  for  the  last  two  approaches  considered  in 


Sec.  IV. 


When  working  with  discrete  probability  laws,  Bigelow  [4]  has 
suggested  an  alternative  approach  that  approximates  a  derivative  by  using 
differences.  The  advantage  of  such  an  approach  is  that  the  normal 
approximation  is  not  needed.  The  disadvantages  are  that  the  resulting 
equations  are  more  cumbersome  and  may  not  assure  more  accuracy  than  simply 
using  t^e  normal  approximation.  Consequently,  the  normal  approximation 


will  be  used  in  what  follows. 
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STATEMENT  OF  CONDITIONS 

Because  the  natural  logarithm  2(»)  is  strictly  monotonica 1 ly 

increasing,  the  objective  function  A  can  be  replaced  with  2(A),  since  a 

solution  maximizes  A  if  and  only  if  it  also  maximizes  2(A).  This 

replacement  will  be  made  in  order  to  convert  the  product  expression  in 

the  objective  function  (2)  into  a  summation.  A  Lagrange  multiplier  v 

will  be  used  to  include  the  budget  constraint  (8)  into  the  modified 

objective  function.  Thus  the  problem  of  maximizing  2(A)  becomes: 

determine  the  stock  levels  x..  to  maximize 

ije 

I  J(i)  2 


STK . .  ] 

i  J  e 


L  =  2(A)  +  vV  Y  y  V.c.  .  [NB  .  .  x.  . 

l  l j  lje  ij( 

i=l  j=0  e=l 

subject  to  the  constraints 


LK. .  <  (x. .  -  m. .  )/S(m. .  )  <  UK . .  . 

lje  lje  lje  lje  lje 


For  the  value  x..  to  be  optimal,  this  stock  level  must  satisfy 
lje  ’  3 

one  of  the  following  four  conditions  necessary  for  optimality: 


x  .  .  =  STK . .  /NB . 

lje  lje  le 


x..  =  S  ( m  .  .  )LK..  +  m ,  . 

ije  ije  l j  e  lje 


x .  .  =  S (m . .  )UK .  .  +  m  .  . 

ije  ije  ije  ije 


3  2  (A)  /  3x  .  .  +  vN'B  .  V.c..  =  0.  (14) 

i J e  le  l  ij 

Constraint  (11)  is  a  lower  bound  because  of  existing  assets,  (12)  is  a  lower 
bound  on  safety  stock,  (13)  is  an  upper  bound  on  safety  stock,  and  (14) 
must  be  satisfied  if  none  of  the  preceding  bounds  are  active. 

It  is  convenient  to  define  the  normalized  stock  level 


k..  =  (x  .  .  -m..)/S(m.,). 

i  je  ije  ije  ije' 


When  using  the  normal  approximation  for  the  pipeline  distribution,  it 
can  be  shown  [6,  p.  446]  that  the  expected  backorders  can  be  computed  as 


EBO. .  =  NB .  S (m . .  )[0(k. .  )  -  k. .  $(k . .  )] , 

lje  le  lje  lje  lje  lje  1 


where  the  normal  density  is  given  by 


0 (r )  =  ( 1/ 2tt)  1/2  exp(-r2/2) 


and  the  right  tail  of  the  distribution  function  by 


SO 


$  (k)  -  J’  0  ( r ) dr . 

k 


In  the  derivations  that  follow,  two  key  results  are  used  repeatedly: 
the  derivative  of  $(r)  is  -0(r),  and  the  derivative  of  0(r)  is  -r0(r). 

Equation  (14)  will  be  written  in  detail  for  each  possible  case. 

First  consider  the  partial  derivative  of  the  Lagrangian  Eq .  (9)  with 
respect  to  the  LRU  stock  level  at  a  base.  According  to  the  chain  rule  for 
derivatives , 


32(A)/3x.01  =  (l/Qi0)OQ.0/aEB0.01)(aEB0.01/3x.01), 


('  1 


where 


3qio/aEBOioi  -  E  (vN)3w3EBOi 

k=l 


iOl 


3qki0/3EB°i01 


(a  b  ,U,  ,/t.  ,n)(l  -  U.  .EBO  nl/t,  .n)akiO_1 
kiO  ki  ki  kiO  ki  iOl  kiO 


and 


3EBOi01/3xi01  =  -NBn#(kioP 


Consequently  Eq .  (14)  becomes 


-Ci/Qio)aQ±0/3EBOio1#^kio1)  +  vVicio  =  °> 


and  ^Q^o^3EB<^i01  ‘*'S  &iven  a^ove  ■ 


Next  consider  the  partial  derivative  of  the  Lagrangian  Eq.  (9)  with 
respect  to  the  LRU  stock  level  at  the  depot.  According  to  the  chain 
rule , 


3fc(A)/3x 


i02 


(1/Qi0)(3Qi0/SEB0i01)OEB0i0l/9n'i01)(3mi01/3xi02)> 


where 


3EBOi01/3mi01  =  NBil[#(ki01^  +  ^ki01)S'(mi01^ 


3mi01/3xi02  =  -$(kiQ2)/NBil 


and  S '  ( • )  is  the  derivative  of  S ( * ) .  Because  NB^  =  1>  Eq.  (14)  becomes 


-U/Q10)SQi0/3EB0io.|f(kioi) 


+  *Ck101)S’(i»i01))*(k102)  +  W.c.0  =  0. 


Next  consider  the  partial  derivative  of  the  Lagrangian  Eq.  (9)  with 
respect  to  the  SRU  stock  level  at  a  base.  According  to  the  chain  rule, 


3H(A)/3x  = 


( 1/Q,’ n )  ( n  / ®EB0  .  _  . )  (3EBO  n  /3m  )  (3m  n  /3x,  . , )  , 


where 


,m101/Sxijl  = 


and  the  other  derivatives  were  given  earlier.  Thus  Eq.  (14)  becomes 


•(1/QiO>3Qio/9EBOioi[^kiOl) 


+  ♦(ki01>s,(®i01))*(kijl>  +  vVicii  =  °- 


(18) 


And,  finally,  consider  the  partial  derivative  of  the  Lagrangian 
Eq.  (9)  with  respect  to  the  SRU  stock  level  at  the  depot.  According  to 
the  chain  rule, 


a*(A)/3x  = 


(l/Q10)(SQi0/SEB0i0i)OEB0i01/3„.0i,(3m.01/Sxlj2). 


where 


S"i01/3xij2  ’  lSEBOijl/axij2  +  (1  -  0IMij)aEB0ij2/axij2l/NBil 


9EB0. . ,/9x. ..  =  ( 9EB0 . . , / 3m . . , ) ( 9m . . ,  / 9x  . . „ ) 
ijl  i  J  2  ijl'  ijl'  ijl'  ij2 


9£B°i  j  i  /  9rn  i  j  i  “  NBii^Ck^i)  +  0(k..1)S,(tn..1)] 


ij- 


ij- 


ij 


a”ijl/axij2  "  -#<kij2)0!Kij/NBi. 


aEBOij2/axij2  -  -»<kiJ2)' 


Consequently  Eq.  (14)  becomes 


-d/QioiaQio/a^ioit^kioi)  +  *(ki01)s’ (m.01)]{(i  '  0INij)f (kij2) 

+  [$  (k .  .)  +  0(k.  )S'(m.  )  ] $(k .  )OIM, .}  +  vV.c..v  =  0. 

ijl  ijl  ijl  ij2  ij  i  ij 
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DIFFERENTIATION  OF  CONDITIONS 

Define  the  column  vectors  X.  =  (x . .  )  and  M.  =  (m . .  ).  Let  Z 

1  ije  1  ije 

represent  an  arbitrary  independent  variable.  Suppose  that  the  AAM  has 
been  solved  for  a  nominal  case,  corresponding  to  the  value  Z^,  thereby 
determining  the  nominal  values  X^  and  M  ^  for  the  two  vectors.  For  the 
purpose  of  developing  ORACLE'S  database,  it  is  desirable  to  compute  the 
partial  derivatives  3X^/3Z  and  3M^/3Z,  as  evaluated  at  the  nominal 
values.  This  subsection  will  show  how  these  derivatives  can  be  computed 
by  differentiating  the  necessary  conditions  for  optimality  and  mean 
pipeline  equations. 

For  each  combination  of  indices  (i,  j,  e) ,  the  previous  subsection 
showed  that  there  is  a  necessary  condition  for  optimality  in  the 


r  ,7o  v  0  M  0  o,  _  . 

F.  .  (Z  ,X .  ,M .  ,v  )  =  0, 
ije  l  i 


where  F..  is  either  (11),  (12),  (13),  or  (14),  and  Eq.  (14)  is 
i  J  e 

either  (16),  (17),  (18),  or  (19).  Equations  (18)— (19)  couple  the 
variables  for  an  SRU  with  those  for  its  associated  LRU. 

Also  for  each  combination  of  indices  (i,  j,  e) ,  Sec.  II  showed 
that  there  is  a  mean  pipeline  equation  in  the  form 


G.  .  (Z° ,X . ° ,M . °)  =  0, 
ije  l  i 


where  G..  is  either  (3),  (4),  (5),  or  (6),  and  the  expected  backorders 
1 J  ® 

in  these  expressions  are  evaluated  with  (15).  Equation  (6)  couples  the 
variables  for  an  LRU  with  those  for  all  associated  SRUs . 

Differentiating  (20)  implies 


J(i)  2 

E  E  t(5rij<i/sxist>(siwsz> + 

s=0  t=l 

+  3F .  .  / az  +  3F .  .  /3v  3v/3Z  =  0, 

i  J  e  i  J  e 


and  differentiating  (21)  implies 

J(i)  2 

EV  [OG..  /  3x .  )  ( 3x .  /  3Z)  +  ( 3G .  .  /3m.  J(3m.  /  3Z )  ] 

1  lje  ist  ist'  lje  ist  ist  1 

s=0  t=l 

+  3G . .  / 3Z  =  0. 
ije 

Because  (6)  and  (IS) - (19)  couple  the  variables  for  an  LRU  with 
those  for  its  associated  SRUs,  differentiating  these  equations  will 
couple  the  partial  derivatives  for  an  LRU  with  those  for  its  associated 
SRUs.  It  therefore  follows  that  the  partial  derivatives  for  a  single 
SRU  cannot  be  determined  separately  from  those  for  other  SRUs  associated 
with  the  same  LRU,  or  from  the  partial  derivatives  for  the  LRU  itself. 

In  other  words,  (22) -(23)  represent  a  system  of  simultaneous  linear 
equat ions . 

If  we  consider  only  LRU  i  and  assume  for  the  moment  that  3v/3Z  is 

known,  then  there  are  2[J(i)+l]  unknown  variables  3x . .  / 3Z  plus 

lje  r 

2[J(i)+l]  unknown  variables  3m..  / 3Z,  or  4[J(i)+l]  unknown  variables 

1JG 

altogether.  Because  (22)-(23)  represent  a  system  of  4[J(i)+l] 
simultaneous  linear  equations,  there  are  sufficient  equations  to  compute 
a  unique  situation. 

Because  F^g  is  a  linear  function  of  the  Lagrange  multiplier  v, 

(22)  will  be  independent  of  the  nominal  value  v°  for  this  multiplier. 
Thus  it  is  not  necessary  to  know  this  nominal  value  to  compute  the 
partial  derivatives  of  the  stock  levels  and  pipelines.  This  latter 
observation  is  important,  because  the  AAM  does  not  use  a  Lagrange 
multiplier  solution  procedure,  so  the  nominal  value  for  this  multiplier 
may  not  be  known. 

Define  the  column  vector  Y.  to  consist  of  both  X.  and  M..  It  is 

i  ii 

convenient  to  rewrite  (22) - (23)  as 


H.3Y./3Z  +  E.  +  R.3v/3Z  =  0, 

ii  ii 


where  each  element  of  the  matrix  H.  is  either  3F . .  /3x.  ,  3F . .  /3m,  , 

l  lje  ist  lje  ist 

3G . .  /3x.  .  ,  or  3G .  .  /3m.  _ :  each  element  of  the  matrix  E.  is  either 
lje  ist.  lje  ist’  i 

3F . .  / 3Z  or  3G . .  / 3Z;  and  each  element  of  the  matrix  R.  is  either 


3F  .  .  /3v  or  zero.  The  solution  to  this  system  of  simultaneous  linear 
lje  J 


equations  can  be  obtained  by  first  inverting  the  matrix  and  then 


computing 


3Y./3Z  =  H.  [  -E  .  -  R,3v/3Z1 . 
1  11  1  J 


The  appropriate  value  for  the  partial  derivative  3v/3Z  can  be 


determined  in  the  following  way.  Define  the  row  vector  such  that  the 


product  C^3Y^/3Z  is  equal  to 


J(i)  2 

y  y  NB.  V  .  c .  .  3x ,  .  / 3Z . 
Z— <  le  l  ij  ije 

j=0  e=l 


If  the  total  budget  B  remains  fixed,  then  the  budget  constraint  (8)  implies 


Y  Ci3Yi/3Z  =  0. 


After  substituting  (25)  into  (26),  we  find  that 


Y  Ci<Hi"1t-Ei  '  R^v/SZ]}  =  0, 


implying  that 


3v/3Z  = 


1  ,  1 

Y  C.H._1E.  /  Y  C.H._1R.  . 

Z— t  ii  l  /  Z—/  ii  i 


By  substituting  (27)  into  (25),  the  desired  partial  derivatives  of  the 


stock  levels  can  be  obtained. 


In  summary,  the  algorithm  for  computing  the  partial  derivatives  of 
the  stock  levels  with  respect  to  an  arbitrary  independent  variable  has 
the  following  steps: 

1.  First  solve  the  AAM  for  the  nominal  case,  thereby  determining 

the  nominal  vector  values  X . ^  and  M.^  for  i  =  1,  ...,  I: 

x  1 

2.  By  differentiating  the  necessary  conditions  for  optimality, 

differentiating  the  mean  pipeline  equations,  and  then 

evaluating  these  relations  at  the  nominal  values,  construct 

the  matrices  H.,  E.,  and  R.  for  i  =  1,  ....  I; 
i  i  l 

3.  Invert  each  matrix  H_^,  which  is  equivalent  to  solving  a 
system  of  linear  equations  with  4[J(i)+l]  unknowns  for 

i  =  1,  • ■ I; 

4.  By  using  (27),  compute  the  partial  derivative  of  the  Lagrange 
multiplier  with  respect  to  the  independent  variable; 

5.  And  by  using  (25),  compute  the  partial  derivatives  for  the 
stock  levels  with  respect  to  the  independent  variable. 

The  main  computational  work  in  this  algorithm  seems  to  be  inverting 
the  matrix  for  each  LRU  i,  which  is  step  3.  However,  this  matrix 
has  a  special  structure  that  may  be  exploited.  In  particular,  this 
matrix  may  be  written  as 


Hll 

H12 

0 

0 

H21 

H22 

H23 

H24 

H31 

H32 

H33 

H34 

0 

0 

0 

H44 

where  column  1  corresponds  to  base  stock  levels,  column  2  to  base 
pipelines,  column  3  to  depot  stock  levels,  and  column  4  to  depot  pipelines 
and  the  rows  are  similarly  defined.  Those  portions  of  the  matrix  that 


have  nonzero  elements  are  indicated  by  the  symbol  HIJ,  where  I  and  J  can 
be  1,  2,  3,  or  4.  Because  H44  is  the  only  nonzero  submatrix  on  the 
fourth  row,  it  follows  that  the  overall  matrix  1L  can  be  inverted  in  two 
stages:  First  invert  H44,  and  then  invert  the  matrix  consisting  of  the 
submatrices  HIJ  for  I  and  J  between  1  and  3.  However,  the  inversion  of 
H44  is  trivial  because  it  is  a  diagonal  matrix.  Thus  the  amount  of 
work  required  to  invert  H^  can  be  reduced  to  solving  a  system  of 
simultaneous  linear  equations  with  only  3[J(i)  +  1]  unknowns. 

It  is  possible  to  exploit  the  structure  of  H^  even  further  by 
developing  an  inversion  algorithm  that  takes  advantage  of  the  zeroes  on 
the  first  row.  However,  such  an  algorithm  would  be  more  complex  than 
the  one  just  described  because  Hll  is  not  diagonal  and  H12  is  not  zero. 

For  simplicity,  Sec.  II  gave  the  equations  for  the  AAM  for  the 
case  in  which  there  is  only  one  level  of  indenture  and  in  which  an  SRU 
can  be  associated  with  only  one  type  of  LRU.  However,  the  new  version 
of  D041  will  allow  several  levels  of  indenture,  allow  an  SRU  to  be 
associated  with  several  types  of  LRUs,  plus  use  a  product  relationship 
(rather  than  a  summation)  for  the  SRU  expected  backorders  in  (6). 

Because  Eqs .  (20)  and  (21)  are  sufficiently  general  to  cover  the 
latter  changes,  the  foregoing  algorithm  is  also  applicable  to  the  new 
version  of  D041. 

Actually  developing  this  algorithm  would  be  quite  time  consuming. 

Such  development  must  include  program  design,  coding,  and  validation.  The 
main  difficulties  would  be  the  complex  matrix  manipulations,  including 
construction  of  the  matrices,  inversion,  and  multiplication. 


IV.  ALTERNATIVE  ORACLE  ESTIMATION  METHODS 


The  indenture  relationships  make  the  approach  of  the  last  section 
difficult  to  implement,  because  these  relationships  result  in  a  set  of 
simultaneous  linear  equations  that  must  be  solved  for  each  LRU  in  the 
weapon  system.  Consequently,  this  section  will  discuss  four  simpler 
alternative  methods.  The  first  method  treats  the  SRUs  as  though  they  were 
LRUs,  eliminating  the  indenture  relationships.  The  second  method  retains 
the  indenture  relationships  but  in  a  simplified  form.  Either  of  these 
approximate  methods  would  be  used  only  to  compute  the  partial  derivatives, 
as  the  nominal  values  would  still  be  obtained  by  solving  the  AAM  without 
any  additional  approximation.  The  third  method  is  based  upon  making 
multiple  runs  of  the  AAM  for  the  same  year.  The  fourth  uses  regression 
analysis  to  analyze  data  obtained  for  several  years. 

ELIMINATION  OF  INDENTURE  RELATIONSHIPS 

For  the  purpose  of  computing  partial  derivatives,  first  consider 
the  approach  of  treating  the  SRUs  as  though  they  were  LRUs.  In  other 
words,  all  items  (indentured  or  not)  are  treated  as  though  they  were 
installed  directly  on  the  aircraft. 

Bigelow  has  also  recommended  this  approximation  scheme  but  suggested 
that  the  data  be  modified  in  the  following  ways:  The  time  spent  working 
on  an  LRU  to  isolate  the  failed  SRUs  should  be  added  to  the  present  SRU 
repair  times;  the  price  of  an  LRU  should  be  reduced  by  the  purchase  price 
of  a  full  complement  of  SRUs;  and  the  SRUs  installed  on  LRUs  should  be 
counted  as  part  of  the  available  SRU  assets. [3]  Bigelow's  suggested 
modifications  have  several  disadvantages:  Additional  data  are  required 
that  may  not  be  available;  the  residual  price  of  the  parent  LRU  may  become 
small  or  even  negative;  the  resulting  necessary  conditions  for  optimality 
for  the  LRUs  are  different  from  those  given  in  Sec.  Ill  and  thus  have 
different  solutions;  and  in  the  absence  of  extensive  cannibalization,  it 
is  unrealistic  to  consider  the  installed  SRU  components  as  being  assets 


that  are  available  for  repairing  other  LRUs.  Consequently,  the  approach 
considered  here  is  to  treat  the  SRUs  as  though  they  were  LRUs,  but  without 
making  any  of  the  foregoing  data  modifications,  which  in  turn  avoids  all 
of  the  foregoing  disadvantages. 

To  derive  the  formulas  for  this  approximation,  it  is  necessary  to 
generalize  the  definitions  given  in  Sec.  II.  Define 

a,  . .  =  quantity  per  application  of  SRU  (i,j)  on  series  k. 

K.1J 

Thus  the  number  of  units  of  SRU  (i,j)  installed  on  aircraft  of  series  k  is 

t,  .  .  =  a.  .  .  b.  .  n.  , 
kij  kij  ki  k 

where  b^  and  n^  were  defined  in  Sec.  II.  The  probability  that  an 
aircraft  of  series  k  is  not  missing  a  unit  of  SRU  (i,  j)  is 

Okij  =  1  •  \i  +  bki(1  -  \iEBOijl'tkij>akiJ, 

where  U^  was  defined  in  Sec.  II.  The  probability  that  a  randomly 
chosen  aircraft  is  not  waiting  for  a  spare  unit  of  SRU  (i,  j)  is 

K 

Qij  =  £  (nk/N)c*kij’ 
k=l 

where  N  was  defined  in  Sec.  II.  And,  finally,  the  probability  that  a 
randomly  chosen  aircraft  is  not  missing  any  of  its  reparable  components 
is 

I  J(i) 

A=nn  v 

i=l  j=0 


which  replaces  Eq.  (2)  as  the  objective  function. 


The  next  step  is  to  modify  the  mean  pipeline  equations  that  were 
given  in  Sec.  II.  For  i  =  1,  ....  I  and  j  =  0,  ...,  J(i),  the  mean  base 
pipeline  is 

mijl  =  [ X  RTS .  .RT  +  X.  (1  -  RTS. ,)0ST. .  +  EBO . . „ ] /NE . , , 

ljl  ij  ijl  ijlv  ij'  ij  ij  2 1 '  il* 

and  the  mean  depot  pipeline  is 

m. .„  =  X .  [ RT . .  ( 1  -  DC. .)  +  DL. .DC. .] , 

ij2  ij2  i J  2  ij'  ij  ijJ’ 

where  the  other  parameters  were  defined  in  Sec.  II. 

Equations  (11)-(14)  in  Sec.  Ill  still  represent  the  necessary 
conditions  for  optimality,  except  that  Eq.  (14)  can  now  be  simplified.  For 
i  =  1,  ...,  I  and  j  =  0,  ...,  J(i),  the  partial  derivative  of  the  Lagrangian  (9) 
with  respect  to  the  stock  level  at  a  base  is 

- ( 1/Q  .  . )  3Q . . /3EB0 . .  #(k . .  )  +  vV.c..  =  0, 
ij  xij'  ijl  ijly  l  ij 

and  the  partial  derivative  with  respect  to  the  stock  level  at  the  depot  is 

-(i/q1.)aQ../3EB0.jl[i(kijl)  +  ♦<k1Jps1(-1J1)|i(k1J2)  *  W.C..  =  0. 


where 


!ij/aEBOiji  -  E  <vN>avySEBoiji 


3^kij/aEBOiji  ■  -(v  b  Ai'bij)'1  ■ 


Let  Z  be  an  arbitrary  independent  variable.  The  approach  of 
Sec.  Ill  can  be  used  to  compute  the  partial  derivatives  of  the  stock 
levels  and  mean  pipelines  with  respect  to  Z,  namely  3x^^/3Z, 

3m.. ,/3Z,  and  3m. ,/3Z.  This  approach  requires  that  a  four-by-four 


system  of  linear  equations  be  solved  for  each  item  (i,  j),  where  these 
equations  art.  obtained  by  differentiating  the  foregoing  necessary 
conditions  for  optimality  and  mean  pipeline  equations.  However,  these 


equations  imply  that  3nUj9/3Z  can  be  evaluated  independently  of  the 


other  variables,  and  that  3m_^/3Z  can  be  expressed  as  a  linear  function 


of  3x;.„/3Z.  Thus  through  substitution,  the  four-by-four  set  of 


equations  for  each  item  can  be  reduced  to  only  a  two-by-two  set,  leaving 


only  the  variables  3x,.,/3Z  and  3x.../3Z  to  be  determined.  In  other 

ijl  ij2' 


words,  treating  SRUs  as  LRUs  can  reduce  the  main  computational 
effort  of  the  Sec.  Ill  algorithm  to  solving  only  a  two-by-two  system 
of  equations  for  each  item,  which  basically  is  the  same  approach  that 
was  described  by  Bigelow  for  minimizing  expected  backorders . [4] 

Of  course,  treating  the  SRUs  as  LRUs  is  an  approximation.  Its 
disadvantage  is  that  there  will  be  some  inaccuracy,  while  its  advantage  is 
that  the  resulting  algorithm  will  be  much  easier  to  program  and  faster  to 
run  compared  with  the  more  exact  approach  considered  in  Sec.  III. 


SIMPLIFICATION  OF  INDENTURE  RELATIONSHIPS 

The  second  approach  toward  simplifying  the  Sec.  Ill  algorithm 
is  to  retain  the  indenture  relationships  but  in  a  simplified  form.  In 
particular,  the  following  approximations  are  made  for  the  purpose  of 
computing  the  partial  derivatives: 


1.  The  normal  distribution  is  used  as  the  probability 
distribution  for  each  pipeline; 

2.  There  is  a  single  level  of  indenture; 

3.  After  setting  OIIL^  =  0  for  each  SRU  (i,  j),  the 
mean  base  pipelines  for  SRUs  and  LRUs  are  given  by 
Eqs.  (5)-(6) ; 

4.  In  the  base  case,  conditions  (11)-(13)  are  not  active 
for  LRU  depot  stock  levels. 

Approximation  1  was  also  made  in  Sec.  Ill  for  the  formulas  for 
computing  the  partial  derivatives  even  to  be  defined. 


Although  approximation  2  was  also  used  in  Secs.  II  and  III,  the 
AAM  does  allow  several  indenture  levels  and  the  algorithm  considered  in 
Sec.  Ill  can  be  extended  to  allow  several  indenture  levels.  But,  and 
this  point  must  be  emphasized,  having  a  single  indenture  level  is 
required  by  the  approach  considered  here. 

The  AAM  already  makes  the  approximation  that  the  calculation  of  depot 
LRU  backorders  can  be  uncoupled  from  that  for  depot  SRU  backorders.  But 
in  an  analogous  way,  approximation  3  also  uncouples  the  calculation  of 
base  SRU  backorders  from  that  for  depot  SRU  backorders  with  the  effect 
that  each  additional  SRU  depot  backorder  causes  an  additional  LRU  depot 
backorder  or  an  additional  SRU  base  backorder,  in  either  case  adding  one 
unit  to  the  mean  LRU  base  pipeline.  Although  this  approximation  will 
overestimate  the  mean  LRU  base  pipeline,  it  may  be  acceptable  for  the 
purpose  of  computing  the  partial  derivatives  of  the  stock  levels. 

Approximation  4  requires  that  the  LRU  depot  stock  levels  for  the 
base  case  exceed  the  lower  bounds  because  of  existing  assets  and  safety 
stocks,  while  being  less  than  the  upper  bounds  because  of  safety  stock. 
However,  these  bounds  are  allowed  to  be  active  for  LRU  base  stocks,  SRU 
depot  stocks,  and  SRU  stocks. 

Equations  (11)-(14)  still  represent  the  necessary  conditions  for 
optimality,  except  that  approximation  3  enables  (14)  to  be  simplified  in 
one  case.  The  partial  derivative  of  the  Lagrangian  (9)  with  respect  to 
the  SRU  stock  level  at  the  depot  is  now 


-(1/Qio)3Qi0/3EBOioi[^kiOl)  +  ^ki0l)S'(mi0l)^(kij2 


)  + 


vV.c.  . 
i  ij 


=  0, 


(28) 


whereas  the  partial  derivatives  for  the  other  cases  are  still  given  by  Eqs . 
( 16 ) ■ ( 18 ) . 

For  each  LRU  i,  approximation  4  implies  that  Eq.  (17)  is  the  necessary 
condition  for  optimality  corresponding  to  the  depot  stock  level,  which 
means  that 


♦(ki02>  -  »VicioW<,‘)iO/,EBOi01<*<ki01> 


*(ki01)S'("i01)l)' 


(29) 


29 
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Suppose  that  the  necessary  conditions  (11)-(13)  are  not  active  for  the 
stock  level  of  SRU  (i,  j)  at  a  base.  Then  Eq.  (18)  must  be  one  of  the 
necessary  conditions  for  optimality,  implying  that 

=  ”VijW(3‘W3EBtW#(ki01)  +  *(ki01)S'(”i01)l>' 

Substituting  Eq.  (29)  implies  that 


After  differentiating  this  expression  with  respect  to  an  arbitrary 
independent  variable  Z,  we  have 


,(kijl!3kijl/aZ  -  (V'iO»*(ki02>:ilW8Z’ 


which  expresses  3k..,/3Z  as  a  linear  function  of  3k.„./3Z. 

ljl  i02 

Next  suppose  that  condition  (11)  is  active  for  the  stock  level  of 

SRU ( i , j )  at  a  base,  implying  that  3xi^1/9Z  =  0.  Because  the  normalized 

stock  level  is  defined  as  k...  =  (x  .  . ,  -  m .  . , ) /S (m .  . , ) ,  it  must  be  true 

ljl  ljl  ljl  ljl 

that 


3k 


ijl/3Z  =  '3mijlZ3Zf1  + 


(30) 


Because  approximation  3  implies  that  3m..,/3Z  can  be  evaluated 

ljl 

directly  as  a  known  function  of  Z  and  the  base  case  solution,  Eq .  (30) 

implies  that  3k..,/3Z  can  also  be  evaluated  directly, 
i  J  1 

Finally,  suppose  that  either  condition  (12)  or  (13)  is  active 

for  SRU  (i,  j)  at  a  base.  In  this  case,  the  normalized  stock  level  k.., 

ijl 

is  equal  to  a  constant,  implying  that  the  partial  derivative  3k^.^/3Z 
equals  zero. 

In  summary,  deperding  upon  which  of  the  conditions  (11)-(14)  is 

active,  3k..,/3Z  can  be  expressed  as  a  linear  function  of  3k.„  /3Z,  as  a 
ijl  i02 

known  function  of  Z,  or  as  zero.  However,  the  last  two  options  are 
special  cases  of  the  first  option,  namely  that  3k^^/3Z  can  be  expressed 
as  a  linear  function  of  3kiQ2/3Z.  By  employing  the  same  argument,  it 


can  be  seen  that  the  partial  derivative  3k.. „/3Z  can  also  be  expressed 

1  j  2 

as  a  linear  function  of  3k^9/3Z. 

It  follows  from  differentiating  (15)  that 


3EB0  .  .  ,/3Z  =  NB.  ,,{-S(m.  ,,)§(k.  ,,)3k.  ,,/3Z 
ljl  ljl  ljl  ijl  ijl' 


+  s-Cm±jl>amiJ1/az[*Ck±jl)  - 


(31 


As  already  noted,  3m,.,/3Z  can  be  evaluated  directly,  and  3k.., /3Z  can 

ijl  ijl 

be  expressed  as  a  linear  function  of  Thus  (31)  implies  that 

3EBO..,/3Z  can  be  expressed  as  a  linear  function  of  3k.-„/3Z.  Using  the 
ijl  i02 

same  argument,  both  SEBO^^/^Z  and  SEBCh^/^Z  can  be  expressed  as  a 

linear  function  of  And,  finally,  (6)  implies  that  3m._^/3Z 

can  also  be  expressed  as  a  linear  function  of  Sk^^/^Z. 

Thus,  for  each  LRU  i,  only  two  partial  derivatives  need  to  be 

determined,  namely  Sk^^/SZ  and  Once  Sk^^/^2  is  determined, 

then  3k..,/3Z  and  3k..„/3Z  can  be  evaluated  using  the  previously 
ijl  ij2 

described  linear  relationships.  Because  3m._„/3Z,  3m..,/3Z,  and 

r  i02  ijl 

3m i j ^ / 3Z  can  all  be  evaluated  directly,  3mi0^/3Z  can  be  computed  by 
differentiating  (6).  After  the  foregoing  derivatives  have  been 
determined,  then  the  partial  derivatives  of  the  stock  levels  can  be 
evaluated  with 


3x .  .  / 3  Z  =  S (m  .  .  ) 3k  .  .  /3Z  +  fl  +  k..  S'(m..  ) ] 3m .  .  / 3Z . 
i  J  e  lje  lje  ije  lje  ije' 


for  i  =  1,  I,  j  =  0,  J(i),  and  e  =  1,  2. 


The  column  vectors  X.  =  (x . .  ),  M.  =  (m . .  ),  and  Y.  =  (X.,  M.)  were 

l  ije  l  ije  l  l  l 

defined  in  Sec.  III.  Also  define  =  (k^^,  k^  )  .  indicated  by 

the  preceding  discussion,  approximations  1-4  enable  (24)  to  be 
simplified  as 


H.  ’ 3K./3Z  +  E . '  +  R.  ' 3v/3Z  =  0 

li  l  l 


3Y./3Z  =  T.'3K./3Z  +  P.'  , 


where  the  matrices  H.',  E.',  R.',  T.',  and  P.'  are  all  functions  of  the 
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nominal  case  values,  and  H/  is  a  two-by-two  matrix.  The  solution  to 
this  system  of  linear  equations  can  be  obtained  by  first  inverting  H.' 
and  then  computing 

3Y_./3Z  =  -T. ’{ (H.' )"1[E.'  +  R.’3v/3Z]}  +  P  (32) 

The  appropriate  value  for  the  partial  derivative  3v/3Z  can  be  determined 
by  substituting  (32)  into  (26) . 

Consequently,  approximations  1-4  enable  the  computational  effort  to 
be  greatly  reduced,  because  the  main  effort  is  solving  only  a  two-by-two 
system  of  linear  equations  for  each  LRU  in  the  weapon  system.  The  earlier 
approximation  method,  which  treated  the  SRUs  as  LRUs,  required  that  a 
two-by-two  system  of  linear  equations  be  solved  for  each  item  in  the 
weapon  system.  However,  both  of  these  approximation  schemes  should  be 
contrasted  with  the  more  exact  algorithm  considered  in  Sec.  Ill,  which 
required  that  a  system  of  linear  equations  having  3[J(i)+l]  variables  be 
solved  for  each  LRU  i  in  the  weapon  system. 

MULTIPLE  FU"!S 

Another  approach  to  constructing  the  ORACLE  database  is  based  upon 
making  multiple  runs  of  the  AAM.  This  approach  is  much  simpler  and  easier 
to  implement  than  the  algorithms  considered  earlier. 

For  a  given  weapon  system,  Sec.  I  defined  Z^  to  be  the  mth 
independent  variable  and  to  be  the  nth  dependent  variable,  for 

m  =  1,  . . . ,  M  and  n  =  1,  . . . ,  N.  Let  Z  ^  be  the  value  of  the  mth  independent 

m 

variable  for  the  nominal  case,  and  let  D  ^  be  the  corresponding  value  of 

the  nth  dependent  variable  as  computed  by  the  AAM.  The  ORACLE  database 

consists  of  the  collection  of  partial  derivatives  3D  / 3 Z  ,  as  evaluated 

n  m 

at  the  nominal  case. 

In  this  alternative  method  of  constructing  the  ORACLE  database  for 
a  given  weapon  system,  the  first  step  is  to  specify  a  nonzero 
perturbation  £  for  the  mth  independent  variable,  m  =  1,  . . . ,  M.  The  second 


step  is  to  make  M  additional  runs  with  the  AAM,  where  the  value  of  the 
kth  independent  variable  for  the  mth  additional  run  is 


+  £ 

m 


for  k  t  m 


k  =  m. 


In  other  words,  for  the  mth  run,  the  mth  independent  variable  is  the 
only  independent  variable  that  is  changed  from  the  nominal  case.  Let 
D^™  be  the  corresponding  value  of  the  nth  dependent  variable  for  the  mth 
additional  run,  as  computed  by  the  AAM.  The  final  step  is  to  approximate 
the  desired  partial  derivatives  with 


3D  /  3Z 
n  m 


(D  m  -  D  °)/£ 
n  n  m 


(33) 


for  m  =  1,  . . . ,  M  and  n  =  1,  . . . ,  N. 
By  definition 


3D  / 3 Z  =  lim 
n  m  _ 


(D 


-  D  )  /£ 
n  m 


m 

if  this  derivative  is  defined.  However,  this  derivative  may  not  be 
defined,  because  of  AAM's  use  of  discrete  probability  distributions  to 
represent  the  pipelines.  If  this  derivative  is  defined,  then  Eq.  (33)  will 
approximate  that  derivative,  provided  that  is  chosen  to  be  suitably 
small.  But  if  this  derivative  is  not  defined,  then  Eq.  (33)  could  still  be 
used,  in  which  case  Eq.  (1)  could  be  interpreted  as  approximating  a 
nonlinear  curve  with  a  secant  rather  than  a  tangent. 

The  above  scheme  requires  M  additional  runs,  one  for  each  independent 
variable.  However,  it  may  be  desirable  to  make  2M  additional  runs  to 
compute  two  estimates  of  the  derivative  for  each  independent  variable,  one 
based  upon  a  positive  perturbation  and  another  based  upon  a  negative 
perturbation.  The  average  of  these  two  estimates  could  then  be  used  in 


the  linear  formula  (1).  Alternatively,  both  estimates  could  be  used  to 
define  a  piecewise  linear  representation,  which  is  likely  to  be  a  more 
accurate  model  than  the  linear  formula  (1)  currently  utilized  for  ORACLE 


REGRESSION  OVER  TIME 

Another  method  for  constructing  the  ORACLE  database  utilizes 
data  obtained  for  several  years.  This  method  is  also  fairly  simple  and 
easy  to  implement. 

Let  be  the  value  of  the  mth  independent  variable  in  year  t,  and 

let  D  ^  be  the  corresponding  value  of  the  nth  dependent  variable  as 
computed  by  the  AAM.  A  linear  formula  is  used  to  relate  the  dependent 
variables  to  the  independent  variables: 

M 

n  t  0  v-**  m_  t 

n  n  L-j  n  m 
m=l 

It  follows  that  a^m  can  be  interpreted  as  being  the  partial  derivative 
of  the  dependent  variable  D  with  respect  to  the  independent  variable 
Z  for  any  year  t. 

The  suggested  approach  is  to  obtain  actual  historical  values  for 
the  relevant  dependent  and  independent  variables  and  then  to  compute 
the  coefficients  a™  by  using  linear  regression  analysis.  However,  such 
an  approach  has  several  difficulties.  First,  the  D041  requirements 
model  has  evolved  over  time,  and  the  AAM  has  not  yet  been  used  during 
any  past  year.  Consequently,  it  would  be  necessary  to  recompute  the 
dependent  variables  for  past  years  with  the  AAM  to  obtain  partial 
derivatives  that  would  be  appropriate  for  the  AAM. 

Second,  the  magnitude  of  a  dependent  variable  in  one  year  may  be 
partly  dependent  on  decisions  made  in  earlier  years.  For  example,  the 
buy  requirement  for  any  year  is  partly  a  function  of  the  percentage  of 
the  previous  year's  requirement  that  was  funded.  Because  some  dependent 
variables  are  partly  functions  of  past  history  and  not  just  of  the 
independent  variables,  there  may  be  a  significant  residual  error  in  a 
regression  approach. 
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And  third,  the  values  of  some  independent  variables  may  be 
correlated.  For  example,  the  age  of  a  weapon  system  and  the  number  of 
flying  hours  per  year  may  both  be  monotonically  increasing  functions 
over  time.  This  phenomenon  is  called  multicol linearity  or 
intercorrelation,  and  it  may  produce  significant  errors  when  using  a 
regression  approach. 


V.  CONCLUSIONS 


Five  methods  have  been  presented  for  developing  the  ORACLE  database 
for  a  given  weapons  system:  first,  differentiation  of  the  necessary 
conditions  for  optimality  and  mean  pipeline  equations  (DIFF);  second, 
elimination  of  indenture  relationships  (ELIM) ;  third,  simplification  of 
indenture  relationships  (SIMP);  fourth,  multiple  runs  of  the  AAM  (MULT); 
and  fifth,  regression  over  time  (REGR) .  Table  1  ranks  these  five 
methods  with  respect  to  simplicity,  accuracy,  and  execution  time. 

Because  MULT  requires  only  that  multiple  runs  be  made  of  an 
existing  model,  it  is  judged  in  Table  1  as  being  the  simplest  method  to 
implement.  REGR  is  also  a  fairly  simple  method,  but  it  requires  some 
additional  work:  To  obtain  partial  derivatives  appropriate  for  the  AAM, 
it  is  necessary  to  recompute  the  dependent  variables  for  past  years  with 
the  AAM  and  then  estimate  a  regression  formula.  However,  neither  of 
these  approaches  requires  much  in  the  way  of  new  software  development. 
The  other  approaches,  which  are  based  upon  differentiating  necessary 
conditions  for  optimality,  require  major  efforts  in  software 
development,  as  well  as  changes  in  software  each  time  that  the  AAM  is 
changed . 

Assuming  that  Eq.  (1)  is  used  to  estimate  the  changes  in  dependent 
variables  as  a  function  of  arbitrary  changes  in  independent  variables, 
which  method  of  computing  the  partial  derivatives  would  yield  the  most 
accurate  results?  One  way  of  comparing  these  five  methods  is  to  think 
of  DIFF  as  using  a  tangent  to  represent  a  nonlinear  curve,  ELIM  and  SIMP 
as  using  approximate  tangents,  MULT  as  using  a  secant,  and  REGR  as  using 
an  approximate  tangent.  If  the  curve  is  convex  or  concave,  then  the 
tangent  would  be  the  better  representation  for  a  small  change  in  an 
independent  variable,  whereas  the  secant  would  be  better  for  a  change  in 
an  independent  variable  that  was  larger  than  the  perturbation  used  to 
define  the  secant.  The  relative  accuracy  of  these  methods  depends  upon 
the  changes  chosen  for  the  independent  variables.  However,  to  compute 
their  tangents,  DIFF,  ELIM,  and  SIMP  used  the  normal  distribution  to 


Table  1 


RANKING  OF  METHODS  FOR  DEVELOPING  THE 
ORACLE  DATABASE 
(1  =  best,  5  =  worst) 


Method 

Simplicity 

Accuracy 

Execution 
T  ime 

DIFF 

5 

3 

3 

ELIM 

3 

5 

2 

SIMP 

4 

4 

1 

MULT 

1 

2 

4 

REGR 

2 

1 

5 

approximate  the  probability  distribution  for  each  pipeline.  According 
to  the  Central  Limit  Theorem  (4),  this  approximation  would  be  good  for 
items  with  fairly  high  failure  rates  but  not  for  items  with  fairly  low 
ones.  It  is  not  clear  how  accurate  these  methods  would  be,  even  for 
small  changes  in  the  independent  variables.  Because  MULT  and  REGR  do 
not  need  to  use  the  normal  approximation,  the  latter  approaches  appear 
to  have  a  slight  advantage  with  respect  to  overall  accuracy. 

If  the  goal  of  ORACLE  is  to  predict  values  of  dependent  variables 
during  a  future  year,  then  REGR  has  an  additional  advantage  with  respect 
to  accuracy:  Historical  trends  can  be  estimated  by  using  the  year  as 
one  of  the  independent  variables,  and  Eq.  (34)  can  be  used  to  project 
these  trends  into  the  future.  Because  REGR  is  the  only  method  that  uses 
data  from  several  earlier  years,  it  is  the  only  method  that  can  capture 
historical  trends.  If  sufficient  historical  data  are  available,  REGR 
may  be  the  most  accurate  model  for  forecasting. 

To  compute  the  partial  derivatives  of  the  dependent  variables  with 
respect  to  a  single  independent  variable,  DIFF  requires  that  a  system  of 
simultaneous  linear  equations  having  3(J(i)  +  1)  variables  be  solved  for 
each  LRU  i  in  the  weapon  system;  ELIM  requires  that  a  system  of  linear 


equations  having  two  variables  be  solved  for  each  item  in  the  weapon 
system;  SIMP  requires  that  a  system  of  linear  equations  having  two 
variables  be  solved  for  each  LRU  in  the  weapon  system;  MULT  requires 
that  one  additional  run  of  the  AAM  be  made,  which  means  solving  a 
nonlinear  programming  problem;  and  REGR  requires  that  one  additional  run 
of  the  AAM  be  made  *or  each  past  year  of  data  to  recompute  the  dependent 
variables,  followed  by  estimating  one  regression  formula. 

The  number  of  unknown  variables  is  the  same  for  the  first  four 
approaches,  namely  four  times  the  number  of  items  in  the  weapon  system 
(corresponding  to  the  stock  level  and  pipeline  at  each  echelon  for  each 
item).  However,  solving  systems  of  linear  equations  should  be  faster 
than  solving  a  nonlinear  programming  problem  having  the  same  number  of 
unknown  variables.  Thus,  in  Table  1,  SIMP  is  judged  as  having  the 
smallest  execution  time,  followed  by  ELIM,  DIFF,  MULT,  and  REGR  in  that 
order . 

The  foregoing  remarks  regarding  the  five  methods  are  summarized  in 
Table  1  with  respect  to  the  three  indices  of  simplicity,  accuracy,  and 
execution  time.  A  method  is  said  to  dominate  some  other  method  if  its 
values  for  all  three  indices  are  smaller  than  those  for  the  other 
method.  As  indicated  in  Table  1,  none  of  the  five  methods  dominates  any 
other  method. 

The  only  method  for  which  there  is  some  computational  experience  is 
the  final  one,  REGR.  Volume  II  of  this  report  discusses  the  performance 
of  that  approach  in  the  context  of  requirements  forecasting. 
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